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Abstract 

Two-photon correlations in the thermal radiation field of a double-star system 
are studied. It is investigated how the differential Hanbury Brown and Twiss (HBT) 
approach can be used in order to determine orbital parameters of a binary star. 
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1 Introduction 


Hanbury Brown and Twiss have discovered the remarkable fact that photons, emitted by 
some thermal light-source, tend to arrive at distant detectors in correlated pairs [H Ej - As 
recognized by Hanbury Brown and Twiss, this HBT effect of photon correlation allows 
to determine the angular size of the thermal light-source. Soon afterwards, momentum- 
correlations among two pions created in heavy-ion collisions have been measured by Gold- 
haber et al. |3j, which are a model-independent possibility in order to obtain information 
about the spatial size and the time-evolution of the hadronic fireball. Some typical two- 
pion correlation functions measured by experiments at CERN-SPS or at AGS-Brookhaven 
are shown in [4] and [5], respectively. Ever since, the HBT approach was used extensively 
for the study of relativistic heavy-ion reactions (6|; for a historical review see [7]. Mean¬ 
while, the method was extended not only to detect the size of the source of emission, but 
also the speed of radial expansion [8], the rotational motion [9] and the turbulence [ID] of 
the emitting source, also by means of the differential HBT method [IT| H21 13]. 

In the original work of Hanbury Brown and Twiss the HBT effect has been exploited 
in order to determine the spatial size of stars. Here, in our investigation we will revisit 
that astronomical problem and apply these new approaches to astrophysical situations 
with the main focus on binary systems. Accordingly, each companion of a double star is 
considered to be an individual thermal source of light and the photons of these stars are 
HBT correlated. Using this assumption, we will examine the possible applicability of the 
differential HBT method to determine some orbital parameter of a binary system: the 
semi-major axis A, orbital period T, and orbital speed of the both components v^,Vb', 
see Fig. [I] 

At the moment being, there are about 10 5 binary systems (resolved, astrometric, 
eclipsing, spectroscopic) are known, according to the Washington Double Star Catalog 
H. However, only for a small part of all these binaries the complete set of all seven orbital 
elements have been determined so far. Consequently, the reason for our investigation is 
twofold: (i) The suggested approach could demonstrate the presence of the HBT effect 
in binary systems, (ii) Such an approach could allow for determining orbital parameter 
independently from standard astrophysical or astrometric approaches. 

The paper is organized as follows: In Section [2] the HBT effect is reviewed for basic 
understanding, both in terms of classical electrodynamics and quantum electrodynamics. 
In Section [3] it is discussed how the advanced HBT approach used in heavy-ion physics 
has to be modified for the case of astrophysics. The emission function as essential element 
of the two-photon correlation function is determined in Section [4] Using this modified 
approach, the two-photon correlation function for a binary system is determined in Section 
[5j In Section [6] the case of one thermal light-source (one star) is considered, while the 
case of two thermal light-sources (binary star) is discussed in Section [Tj A summary and 
outlook is given in Section [8j 

Throughout the article, we use Heaviside-Lorentz units: eo — fJ>o — c — h — 1, and 
for the Boltzmann constant ks = 1. Furthermore, the astronomical unit is denoted by 
1 a.u. = 1.49 x 10 11 m. 
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Figure 1: The orbit of a binary system is defined by seven orbital elements: distance 
L between the center-of-mass of the binary star and the observer, semi-major axis x , 
inclination 0 < i < tt, eccentricity 0 < e < 1, eccentric anomaly 0 < E < 2n, periapsis 
0 < uj < 7T and mass ratio M 1 /M 2 . Other familiar elements (e.g. orbital period T or 
orbital speed V\,V 2 ) can be deduced from these seven orbital elements. The elements x, 
e, and M^/Mb uniquely determine the shape of the both ellipses, while i and u determine 
the orientation of the orbital plane (here the [x, y]-plane), and a parameter E determines 
the actual position of both components of the binary, 00 and G. In case of a symmetric 
system where M = Mi = M 2 and r 0 = r 1 = — t* 2 , an idealized, circular orbit may be 
possible. 

2 HBT effect 

In 1954, Hanbury Brown and Twiss na have realized the idea of an intensity interferome¬ 
ter for classical electromagnetic radiation in the radio-band. Such intensity interferometry 
compares intensities of the incoming radiation field rather than amplitudes at different 
detection points. Originally, the basic idea was that intensity-measurements are able to 
improve the angular resolution of radiation-sources and do not account for the phase of 
the electromagnetic fields, hence are insensitive to phase shifts caused by atmospheric 
scintillations. They have first tested their intensity interferometer by determining the 
known angular size of the Sun, and later they succeeded in determining the angular size 
of the radio sources Cassiopaia A and Cygnus A [T5] . 

In a first laboratory desk experiment, Hanbury Brown and Twiss have extended their 
new technology into the optical range of electromagnetic spectrum. They have used a 
light beam from a mercury vapor lamp, which is actually a thermal light source, and were 
able to detect second-order correlations among photons [T], being the discovery of a new 
physical effect. Afterwards, Hanbury Brown and Twiss succeeded with their technique 
in determining the angular size of the star Sirius by detecting photon correlations in the 
optical region [2]. 

While the correlation of radio waves is an effect which can fully be explained by 
classical electrodynamics, the correlation of photons in the HBT experiment initiated a 
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heated debate about the concept of the photon at that time. Later, Glauber succeeded 
with a comprehensive theoretical explanation of the HBT effect in a series of articles 
[M1CE71 HE]. Especially, Glauber (i) introduced the concept of higher-coherence in quantum 
electrodynamics [Mj, (h) determined the explicit quantumfield-theoretical expression for 
coherent states (Glauber states) [IT], (iii) determined the general expression for the density 
operator for coherent and incoherent states [18], and (iv) has shown that the HBT effect 
does not happen for Laser light sources (described by Glauber states), but occurs for 
extended thermal light sources (described by incoherent states) [T 8 J - These articles were 
the birth of a new branch of science, the Quantum optics, for which he awarded the 
Physics Nobel Prize in 2005. 

In general, in HBT measurements one records photons of a thermal light source by two 
detectors: detector A is located at xa and detects incomming photons at time tA, while 
detector B is located at x B and detects incomming photons at time ts- The spatial dis¬ 
tance between both detectors is cLab = \ x a ~ x b\- Depending on the concrete experiment 
under consideration, two specific cases have to be distinguished: correlation measurement 
in space {tA = t b and xa ^ xb) and correlation measurement in time {tA ^ ts and 
xa = Xb)- In our investigation, correlation measurements in space are relevant, that 
means the incoming photons are recorded at the same time, tA = ts, while the distance 
cIab between both detectors is varied. In this way one can determine second-order cor¬ 
relations (HBT effect) of incoming photons, which are used to determine some physical 
parameters of the thermal light-source. Fundamental aspects of the HBT effect will be 
demonstrated by means of an elementary model for the thermal radiation. 

2.1 HBT effect in terms of classical electrodynamics 

In order to describe the HBT effect it is useful to subdivide the surface of a star with 
stellar radius R into pointlike regions, each of which is considered as an emitter of spherical 
waves of thermal light, i.e. pointlike sources of black-body radiation. The entire radiation 
field of the star would finally been obtained by a summation over all of these pointlike 
regions over the whole surface of the star. 

For the description of spherical waves it is advantageous to introduce spherical coordi¬ 
nates x = ( x,y,z ) —* r = {r } 6,(p), i.e. x = r cos 0 cos 0, y = r cos 6 sin 0, z — r cos 6. 
Here the origin of coordinate system is, first of all, assumed to be located at the center of 
a pointlike radiation source, hence r = \x\ being the distance between the pointlike region 
of emission and some point with spatial coordinate x r (later it will be the detector’s 
position). In general, the radiative part of electric field at sufficiently far distances from 
a localized radiation-source, defined by the wave-zone k r 1 (wave-number k is related 
to wavelength A in virtue of k = 2 7 t/A), is given by the following multipole-expansion for 
one individual frequency 04 = k [19]: 

gi(kr—LJkt) 00 l 

E ( r ,*) = i — rz —E E «im 9l m (0 , <f>) + C.C. , (1) 

K r 1=1 m=-l 

where we have assumed that the radiation source has no magnetic moments; c.c. stands 
for complex-conjugate and l > 1 because there is no monopole radiation. Here, g lm = 
(—i) 1 n x Xi m with the unit direction n = r/r and the vector spherical harmonics Xi m = 
l Yi m /yjl (/ + !), the orbital angular-momentum operator is l = r x V, and the spherical 
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harmonics are Y] m . The eigenvalues of the spherical harmonics are l 2 Y] m — 1(1 + 1) Yi m 
and l z Yi m = mYi m . The electric multipole-coefficients are complex-valued numbers, 
aim = \0Lim\ with phase 

In order to get an idea about the magnitude of wave-number k, we recall that a star 
is an almost ideal black body radiator, that means the wavelength of maximal intensity 
of a star is determined by Wien’s law of displacement, which states that the wavelength 
A mav of maximal intensity of a star with surface temperature T is given by A max = 2.9 x 
1CT 3 meter Kelvin/T. The corresponding wave-vector k max = 2 7r/A max . For instance, 
Sirius has a surface temperature of about 9900 K and we obtain A max = 2.92 x 10 -7 meter 
and fc max = 2.15 x 10' meter -1 . Therefore, it would be meaningful to determine the HBT 
effect with photo-detectors which are sensitive in the vicinity around fc max . 



Figure 2: Two pointlike regions at x\ and X 2 of the extended light source and two detectors 
located at xa and xb ■ The vectors are defined by r i2 = X\ — x 2 , cLab = xa — xb and 
x A + x B d A B , T d AB 

To simplify the notation, we assume the photo-detectors to be sensitive only for incom¬ 
ing light rays which are polarized in an arbitrary but fixed unit-direction e, thus we con¬ 
sider the vector-component of the electric field in this direction: £ (r,t) = e • E (r, t) [IE], 
and we introduce the notation g lm = e ■ g lm . Actually, the radiation distribution caused 
by a superposition of incoherent multipoles in the pointlike radiation-source is isotropic 
[19], hence independent of m. However, in HBT experiments one measures individual 
photons, which in the classical electrodynamics correspond to one single radiation-mode 

lr The electric multipole-coefficients cp m used in our article are related to the electric multipole- 
coefficients ole (l, in) used in [19] via ai m = Zq cie (l, m), where Zq is the wave-impedance in vacuum. 
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of the electric radiation-field in ([Tj) , given by: 

Q i(kr-u} k t) 

£lm (t* > — i T O^lm Qlm (fij 4*) "F C.C. . (2) 

rC V 

For HBT effect, it is sufficient to consider the most simple case of two specific spherical- 
wave modes: one mode l\mi originates from the pointlike region at X\ of the star’s surface, 
while another mode I 2 ITI 2 originates from the pointlike region at X 2 of the star’s surface. 
The collective electric held at detector position x A and at time t is then given by a 
superposition of these waves, cf. Fig. [2j 


£ (x 1 ,x 2 ,x A ,t) ^limi (x 1 ,x A ,t) + 

£hrri2 (x 2 ,X A ,t) . 


(3) 


Using the relation for a spherical wave-mode originating at aq towards detector at x A 

X A — X\ 

with momentum k\ A = k\- --, and another spherical wave-mode originating at X 2 

\X A -Xi\ 

X A — X 2 

towards detector at x A with momentum k 2A = k 2 —--, the scalar product takes the 


form 


XA - X 2 \ 


ki A -(x A -x 1 ) = ki \x A - *i| , 
k 2A -{x A -x 2 ) = k 2 \x A - x 2 \ , 


(4) 

(5) 


where the two electromagnetic spherical-waves, originating from region Xi and x 2) are 
given by 


£lim\ (*® 1 , XAjt'j i \oi[ 1 


m 1 


Sh 


m 1 


£l-2m2 {X2,X A , f) i |cr/ 2m2 | Ql 2 m2 


exp [i (ki 

\x A - 

- x x \ 

— LUk-t t + Pki)\ . 

exp [i(k 2 

k\ 

\x A - 

\x A - 
- X 2 \ 

- *l| 

— U>k 2 t + Pk 2 )] 

k 2 

\x A - 

- X 2 \ 



( 6 ) 

(7) 


where the phases are written explicitly. The electric held at detector B is obtained in a 
very same way, where x A is simply replaced by Xb- The correlation function of second- 
order, for these two spherical-waves is dehned by 


C(x i,x 2 ,x A ,x B ) = 


(I (x 1 ,x 2 , X A: t ) I (a?i, x 2 , x B , t)) 

(I (aii, x 2 , x A , t)) (/ (* 1 , ® 2 , x B , t)) ’ 


( 8 ) 


where the intensities I (x A ,t) and / ( XB,t ) are dehned as square of the absolute value of 
electric held and are measured by two intensity detectors located at x A and xb , 


I (x 1 ,x 2 ,x A ,t) = £* (xi,x 2 ,x A ,t) £ (xi,x 2 ,x A ,t) , 
I(xi,x 2 ,x B ,t) = £* (x 1 ,x 2 ,x B ,t) £(xi,x 2 ,x B ,t). 


(9) 

( 10 ) 


Thermal radiation helds are characterized by the fact that the amplitudes cq m huctuate 
independently of one another, that means their phases (fik are randomly distributed as well 
as their absolute values |cq m |. Accordingly, the thermodynamical average of a function 
O in (|8]) implies, hrst of all, an average over the phases and afterwards an average over 
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the absolute values of the amplitudes. The average over the phases resembles the thermal 
average procedure originally introduced by the Einstein-Hopf model |20j . and is given by: 


(O) = 


1 

27T 


2tt 

J ° (+) d( f ■ 

0 


( 11 ) 


By using (J3]) , (|9]) and (10), and performing an average according to ( |TIj ) one obtains for 
the correlation function of second-order for two spherical-waves with equal frequencies 
Wfci = the following expression: 


C(x \,x 2 ,x A ,x B ) 


= 1 + 2 


cq 


imi 


I" Km/ 


(K mp | + \oil 2 m2 1 


COS 


k (|*i - x A \ + \x 2 - x B \ - \x 2 - x A \ - |a;i - x B \) 

( 12 ) 


As stated above, we also have to perform a thermal average over the absolute values of 
the amplitudes in (12), and one easily obtains: the relations: ( ar, mi | 2 )th = (|cq 2 m 2 1 2 )tn 
and (|cfy imi | 2 |cty 2m2 | ) t h = 2 (|cq imi | 2 ) t n • Then, the correlator in (12) hnally yields 


C (x 1 ,x 2 ,x A ,x B ) = 1 + cos 


k (|£Ci - x A \ + \x 2 - x B \ - |*2 - x A \ - |*i - x B \) 


(13) 


There are two important limits for the correlation function in (13), namely the case of 
stars and the case of heavy-ion collisions, cf. Refs, 
consider in more detail. 


. Both limiting cases we shall 


2.1.1 Correlation function in the limit of stars 

In case of a star we have the following limits: L r\ 2 d AB , where L is the distance 
between star and detectors, ri 2 is the distance among two pointlike regions on the star’s 
surface, and d AB is the distance between both detectors, see Fig|2j In these limits, the 
argument of the cosine-function in (13) can considerably been simplified and we obtain 
for stars (for more details see Section |5|): 


Cst ax (x 1 ,x 2 ,x A ,x B ) = 1 + cos ( k 


d AB ■ ri 2 


(14) 


where k = 2 tt/X is the wave-vector of the radiation emitted by the pointlike sources, 
d AB = x A — x B is the vector from detector A to detector B and r\ 2 = *i — * 2 is the 
vector from source-point X\ to the source-point x 2 on the star’s surface. 


The square of absolute value of amplitude is proportional to the energy: | cq r 


Ek ■ 


2 (°° dE k E k e~ E k/ T 

Accordingly, the thermal average is determined as follows: (|cq m | )th = roe -= T and 


/I 1 + Jo dE k E k e k I 2 0 n |2\2 

{\0ilm\ )th — r°o e -E k / T — 2T — 2(|aj m | } th . 
Jo 


JJ° dE k e~ E k/T 




















So far we have considered two pointlike thermal light-sources of the star. The second- 
order correlation function of classical thermal radiation emitted by the entire surface of a 
star is written as follows: 


C(x A ,x B ) 


(■I (x A ,t) I (x B ,t)) 

(-I (x A ,t)) (I (x B ,t)) ' 


(15) 


The opacity of a star (impenetrability of electromagnetic radiation or visible light) is very 
high so that the radiation held originates from the star’s surface only, more accurately the 
light is emitted by the photosphere of the star. Accordingly, we can obtain the correlation 
function in (15) by a summation of the expression (14) over all possible configurations of 
regions aq and x 2 of the entire surface of a star. Consequently, we have to integrate over 
the surface of a star facing the observer (A star = tt R 2 ) as follows: 


C*Star {p^Ai *®b) 



/ d 2 x 2 C (x 1 ,x 2 ,x A ,x B ) 

J -^star 


4 L 2 2 (kRd AB \ 

>?<P AB IP ‘( t j 


(16) 

(17) 


where J\ is the Bessel function of first kind. The result in (JTt]) reflects the Siegert relation 
for thermal light in classical electrodynamics [24], which relates second-order and first- 
oder correlation function. An example for the correlation function in ( |l7| ) is given in 
Fig. [3] for Sirius, repeating in this way the intensity interferometry analysis of Hanbury 
Brown and Twiss (2j. 

In HBT experiments one determines the angular size of the star, 0 = 2 R/L, as seen 
by an observer at distance L, so that the correlator in (JTt]) reads 


^Star (*^A; *®b) 


1 , 16 t2 ( kd AB Q \ 

+ k 2 d? AB e 2 A 2 ) 


(18) 


which agrees, for instance, with the correlation-function given by Eq. (2) in [25]; recall 
k = 2tt/X. The correlator (18) varies as a function of the telescope separation d AB . 
Thus, by varying the separation of the detectors d AB , one may deduce the apparent angle 
0 of a remote star, even if the source is optically not resolvable. In real experiments 
one compares the correlator determined by real measurements with the correlator in (18) 
(fitting procedure). 


2.1.2 Correlation function in the limit of heavy-ion collisions 


In case of heavy-ion collisions (HIC) we have the following limits: L d AB r 12 , where 
L is the distance between fireball and detectors, r 12 is the distance among two pointlike 
regions on the fireball’s surface and d AB is the distance between both detectors, see Figj2} 
In these limits, the cosine-function in (13) can be considerably simplified and is given by 


Cmc{x 1 ,x 2 ,x A ,x B ) = 1 + cos [k 


dAB • r 12 


(19) 


For comparison with the heavy-ion collision results we recall the practical relation q = 
Qab = kd AB /L, because there the distance of the detectors, d AB , is not measured directly. 
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Figure 3: The normalized correlation function Cgtar as given by Eq. (JTt]) for the following 
parameter (Sirius): stellar radius R = 1.2 x 10 9 m, distance between star and detectors 


L = 8.61ightyear = 8.136 x 10 16 m, and wave-vector k = 2.15 x 10'm 
with (18) we deduce an opening angle 0 = 2 R/L = 2.9 x 10~ 8 rad. 


-l 


By comparing 


We recognize that the correlation function in the limits of heavy-ion collisions in (19) 


agrees with the correlation function in the limit for stars (14). Also in case of heavy- 


ion collisions one has to sum over all individual pointlike regions over the entire fireball 
created in the heavy-ion collision process. However, the description of HBT in heavy-ion 
collisions is more complicated than for stars, because the fireball expands rapidly in time 


and space, and therefore a simple two-dimensional integration procedure like in (16) is 
inapplicable, instead a four-dimensional description is necessary. 


2.2 HBT effect in terms of quantum electrodynamics 

The surprising discovery of quantum correlations of second-order in ordinary light by 
Hanbury Brown and Twiss came to pass when they applied photodetectors in the do¬ 
main of the optical spectrum DEI' The theoretical description of a localized absorption 
process of a photon by photodetectors (e.g. photodiodes, photomultipliers, CCD’s) im¬ 
plies the electromagnetic held to be considered as made of by photons instead of classi¬ 
cal waves. One aspect thereof is the symmetrization of the photon wavefunction under 
coordinate-exchange of two indistinguishable photons, a term which, however, does not 
occur in classical electrodynamics. Accordingly, the use of a photon-counter implies that 
the electric radiation-held E ( r,t ) in (jl| must have to be described as a second-quantized 
quantum-held, indicated by a hat: E ( r,t ) |26, 127] : 

= E (+ \r,t) + E^\r : t), (20) 

„ c_,\ e +i(kr-uj k t) oo l 

E ( r,t ) = +i -—-51 51 azm3z m (M), (21) 

^ ' 7 _ 1 7 
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£ t_) (r,t) = 


^—i(k r—ujfc t) oo 


— I 


kr 


i 

^ 9lr, 

1=1 m=—l 


£ 






( 22 ) 


By comparing the second-quantized electric field operator in (20) with the classical electric 


field in Q, one recognizes that there is formally a replacement of the classical amplitudes 
by operators 0 cq m —» a* m and a* m —>• aj m , where aj m is the Hermitian adjoint of a/ m . 
These photon-operators act in the Fock-space of photon states and obey commutator 
relations in angular-momentum space, given by 


i 


mi , a l 2 m 2 


— Uk Si x i 2 8. 


m\m2 5 


[ai imi j 8l 2m2 ]_ 0 , 


t at 

u hmi ) u 'Z 2 m 2 


= 0. (23) 


The operator aj imi creates one photon with angular-momentum l\mi, while the operator 
diimi annihilates one photon with angular-momentum Ipi. In particular, the vacuum 
state is defined by |vac) = 0, and \l\mi) = a\ |vac) represents a one-photon state, 

while |Zimi ^ 2 ^ 2 ) = r «l 2 m 2 l vac ) represents a two-photon state in the angular- 

v 2 

momentum space. Accordingly, the positive-energy part (21) annihilates one photon 


localized at r, 0,0, t in coordinate-space, while the negative-energy part (22) creates one 
photon localized at r, 0, 0, t in coordinate-space. 

Like in the classical case, we assume the detectors are sensitive only to those photons 
which are polarized in an arbitrary but fixed unit-direction e, and consider the vector- 
component of the electric field-operator in this direction: £ (r,t) = e ■ E ( r,t ) |16j . One 
individual mode of the electric field-operator is then given by 


£(r,t) = £ t+) (r,t)+£ [ 1 (r,t) , 


(24) 


^ e +i(kr-u k t) 

^ (r ? + i ~ ttlm fjlm (^5 0) ) 


£ { = 


kr 

3 - i{kr-u k t ) 


kr 


-a, m g‘ lm («, <t>) 


(25) 


(26) 


The field-operator in Eq. (24) is the quantum-field analogon of the classical electric field 
in Eq. ([2]). 

The theory of HBT effect in terms of quantum electrodynamics has been worked out 
by Glauber in a series of articles [sunns]. We shall recall some parts of this theory 
which are relevant for our investigation. 

The initial state before a one-photon detection at time t is denoted by |i), and the 
final state after detection is denoted by |/) which contains one photon less because it 
has been absorbed by the detector located at xa- The final states form a complete set 
of states, £/ [/) (/| = 1. The corresponding amplitude for a one-photon absorption is 
given by the following matrix element: (/|£d + ) (aj A ,t) |i). The probability of this process 
is given by the square of the absolute value of this amplitude (Fermi’s golden rule), 


= 


(/|£( + ' ) (xA,t) |z) . If one wants to determine the total probability that one 


3 A more accurate second-quantization procedure is described in detail in Haem and takes into account 
the spin of the photons (s 2 = —1,0, +1) by introducing the total momentum-operator j = l + s, where 
s is an operator acting in the Hilbert-space of spin-states of photons. However, the exlicit notation of a 
spin quantum-number would not change the basic arguments given. 
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photon is absorbed at time t by an ideal photon detector located at x A and irrespective 
of the final state, then one has to sum over all final states and obtains, cf. Eq. (2.15) in 


{x A ,t) = J2 |(/l^ (+) ( x A,t) \i) * = {i\£ ( } (x A ,t)£ M (x A ,t) \i ), 

/ 

where in the last term a summation over the complete set of final states has been per¬ 
formed. If the initial state |z) is not a pure state but a statistical mixture of states (thermal 
bath of photons), then the initial state is described by a statistical operator p and the 
average in (27) has to be determined by the trace over the density operator: (i\0\i) —* 
Tr (p O). Similarly, in case of two-photon detection at time t by two detectors located 
at x A and xb, the amplitude of this process reads (/|<E + ) ( x B ,t ) (x A ,t) |z), and the 


(27) 


{/|£ t+) (xa, t) £ (+> (x B ,t) | i) 


(o) 

probability of this process is (Fermi’s golden rule) W>f = 

The total probability, that two photons are absorbed at time t by two ideal photon d< 
tectors located at x A and xb, and irrespective of the final state is given by, cf. Eq. (2.17) 
in [16]: 

W t (2) (x A ,x B ) = J2 |(/l^ {+) ( x B,t) £ {+) (x A ,t) \i) 
f 

= (i\£ { ~ ] {x A , t ) (x B , t ) £ {+) (x B ,t) £ i+) (x A , t ) | i) , (28) 

where in the last line a summation over the complete set of final states has been performed. 


The expression in (28) describes basically the correlation of a two-photon absorption pro¬ 


cess, namely one photon is absorbed at detector position x A and at the same instant of 
time another photon is absorbed at detector position x B . In this sense one can call this 
expression correlation function. For a theoretical description of real correlation experi¬ 
ments it is convenient to consider, Erst of all, two pintlike regions located at aq and x 2 
on the star’s surface, each of which emits one spherical photon-mode. Then, we arrive at 
the following correlation function, 

E ( f\£ {+) (*i, * 2 , x B ) £ (+) (aq, x 2 , x A ) \i) 

C(x i,x 2 ,x A ,x B ) = 


f 


E (/1^ (+) (aq, x 2 , x B ) |z) E ( f\£ {+) (aq, x 2 , x A ) \i) 
f f 


2 > 


where we also have introduced a normalized form of (28), cf. Eq. (4.3) in 


(29) 


The 


correlation function (29) determines the rate of coincidences in the photon count rate of 


two detectors at x A and x B , and is the quantum-field theoretical analog of the classical 
correlator ® for two pointlike regions of thermal radiation. In order to determine the 


correlation function in (29), the electric field-operators in (29) can been obtained by a 


second-quantization procedure of the classical expressions in ([3]) - (J7]) with the aid of 
formal replacements of the classical amplitudes by operators cq m —* d/ m and a* lm —* aj m , 

£ (+) (x u x 2 ,x A ,t) = £ (+) (x u x A ,t) +£ {+) (x 2 ,x A ,t ), (30) 

where these two positive-energy field-operators 

«?(+)/ .x exp[+i(/q \x A -Xi\ -u kl t)] 

M \^A — *^11 

S(+) f .x exp [+i (k 2 \x A — x 2 \ — ui k2 1)] 

£ ^12^12 912^2 7 I I * (32) 
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The field-operator at detector positon xb can be obtained by a replacement of xa by xb 


in Eqs. (30) 


(32). In the initial state in the nominator and denominator in (29) there are 
1 


two photons: |i) = l 2 m 2 ). In the nominator in 

v- 


state (/| = (vac|, while in the denominator in (29) the final state is a one-photon state, 


29) the final state is the vacuum 


either (/| = (l\mi\ or (/| = (Z 2 rrz 2 1, depending on which photon has not been detected. 
We assume that both emitted photons have the same energy k\ = k 2 but, for simplicity, 


here we take either l\ ^ l 2 or mi ^ m 2 . Then, by inserting the field-operators (30) into 


(29) we obtain the following expression for the correlation function: 

2 


C (x 1 ,x 2 ,x A ,x B ) = I'h 

where the photon wavefunction [j] is given by 
^12(*1,*2,*A,*b) = 


12 


1 

72 L 


gik\x!-XA\ gifc \x2—xg\ _|_ gjfc \x2—xa\ e ik\x!-XB\ 


(33) 


(34) 


This wavefunction is valid both for heavy-ion collisions and stars, and symmetric under 
exchange of source-coordinates aq x 2 and detector-coordinates Xa -h- x b ■ After an 
experimental measurement the wavefunction collapses and disappears completely. Let us 
also recall, that \k 12 is not an observable, while |'k 12 | can be determined by experiments. 


For the absolute square of photon wavefunction (34) we obtain 


IT 


12 


= 1 + COS 


*( 


Xi - X A \ + \x 2 - X B \ - \X 2 - X A \ - \Xi - X B \ 


(35) 


We recognize, that the absolute square of wavefunction in (35) equals the classical corre¬ 


lation function given by (13). 


The symmetrization of wavefunction in (34) is a fundamental law of quantum me¬ 


chanics and is independent of any specific performance of an experiment. In other words, 
the photons, during their propagation from the star (or fireball of a heavy-ion collision 
process) to the detectors, do not take care about the existence of some detectors any¬ 
where in the universe. For instance, in real HBT experiments in heavy-ion collisions one 
might want to install about 3000 detectors around the fireball in order to minimize the 
loss of photons which have been emitted by the fireball and to improve the statistics and 
accuracy. Finally, one is searching in the experimental data for all two-photon corre¬ 
lations (HBT effect) measured by any two detectors among these 3000 detectors. But 
the photon (or pion) wavefunction is symmetric in the coordinates of the fireball as well 
as in the coordinates of the final two detectors. The same symmetry exists in case of 
HBT measurements for stars. However, different experiments have different geometrical 
configurations. Depending on these specific experimental configurations, see Fig. [4| some 
parts of the wavefunction in (34) become important, while other parts in the wavefunction 


become negligible. For instance, in case of stars the wave-function is given by Eq. (1) in 
pa. while in case of heavy-ion collisions the wave-function is given by Eq. (8) in [22]. 

Actually, more important is the absolute square of the wave-function. Especially, like 
in the classical theory, there are two relevant limits for the absolute value in Eq. (35), 
stars and heavy-ion collisions, which will be considered in the following. 


4 Let us briefly make some comments about this notion. According to a strict mathematical statement 
of Landau and Peierls [28| . there is actually no photon wavefunction in quantum-field theory. On the 
other side, interference effects in double-slit experiments have shown that the probability-density for 
photons have some kind of analogy with electrons, e.g. [29], and have allowed to reintroduce the concept 
of a photon wavefunction in an appropriate meaning, for a review we refer to 
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Figure 4: The extreme limits for the momenta. Left pannel: photon-momenta in case of 


stars r\2 3> cLab- In this limit the wavefunction T 12 in (34) simplifies to in (36). 


Right pannel: photon-momenta in case of heavy-ion collisions r 12 cLab ■ In this limit 

the wavefunction T 12 in (34) simplifies to Tf 2 IC in (39). 


2.2.1 Wavefunction for stars 

In case of a star we have the following limits; see left pannel of Fig. [4} L ^$> r± 2 cIab, 
where L is the distance between star and detectors, ri 2 is the distance among two pointlike 
regions on the star’s surface, and dAB is the distance between both detectors, see Figj2} 
In these limits, photon wavefunction in (34) can considerably be simplified and is given 
by; cf. j221: 


,t. Star __ 

12 “ x/2 L 


Jlk 1 -(xA-x 1 )+k^-(xB-x 2 )] _|_ Q i[k 1 -(x B -x 1 )+k2-(xA-x 2 )] 


(36) 


where the wave-vector Aq is directed from spatial coordinate aq of the star’s surface to 
detector A or B, and wave-vector k 2 is directed from spatial coordinate x 2 of the star’s 


surface to detector A or B. The derivation of Eq. (36) from Eq. (34) can be performed 


like in Section [5) The square of absolute value of this wavefunction is given by 


. 7 , Star 
1 12 


= 1 + COS k 


d A B ■ r 12' 


(37) 


where dAB = xb~ xa is the vector from detector A towards detector B, and Tq 2 = x 2 — aq 

is the vector from surface-point x\ towards the surface-point x 2 of the star, k — |fci| = 

ri2 

|fc 2 |, anc l we have used the relation k 2 — ki = k - which is valid for the geometry of 


stars; see left pannel in Fig. |4| Obviously, (J3T|) equals the classical expression obtained in 

P|). 
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Like in the classical case, see Eq. (16), in order to obtain the correlation function for 


the whole radiation held of a star, one has to sum over all individual pointlike regions 
over the entire surface of the star facing the observer: 


C'star 3Cb) — 


-p— f d 2 x i f d 2 x 2 .' 

■'T-ct.nr *^-^star ^-^-star 


= 1 + 


1 star J A atar 

4 L 2 


Star 

12 


k 2 d\ B R 2 


Jl 


k R dAB \ 


I > 


(38) 


which is obviously the same expression for the correlation function for stars in Eq. (17), 
previously obtained in the treatment in terms of classical electrodynamics. As specific 


example, Eq. (38) reflects the analogon of Siegert relation in quantum electrodynamics, 
first obtained by Glauber, cf. Eq. (10.26) in (18j, and in general valid for incoherent light. 


2.2.2 Wavefunction for heavy-ion collisions 

In case of heavy-ion collisions we have the following limits; see right pannel of Fig. |4} 
L dAB ^ r 12 ) where L is the distance between fireball and detectors, r i2 is the distance 
among two pointlike regions on the hreball’s surface and dAB is the distance between both 
detectors, see Figj2} In these limits, the photon wavefunction in (34) can considerably be 
simplified and is given by; cf. m 


®HIC _ - 

1 12 — ^ 


1 

V 2 


e i[k A -{x A -x 1 )+k B {x B -x 2 )] _|_ e i[k A (x A -x 2 )+k B -(x B ~x 1 )] 


(39) 


where the wave-vector is directed from the fireball's surface towards detector A, the 
wave-vector k B is directed from the fireball's surface towards detector B. The square of 
absolute value of this wavefunction is given by 




HIC 

12 


= 1 + cos k 


dAB ■ r 12 


(40) 


where r\ 2 = x 2 — X\ is the vector from one surface-point X\ of the fireball of HIC to 
another surface-point x 2 of the fireball, k = \Iza\ = |&a| is the absolute value of wave- 


vector of the photons, and we have used k B — kA = k 
of heavy-ion collisions; see right pannel in Fig. |4} 


d 


■AB 


which is valid for the geometry 


We note that (40) equals the corresponding classical expression obtained in (19). Like 


in the classical case, the absolute square of wavefunction in the limits of heavy-ion col¬ 


lisions in (40) agrees with the the absolute square of wavefunction in the limit for stars 


(37). Let us recall here that the fireball created by a heavy-ion collision expands rapidly 


in time and space, hence the two-dimensional integration procedure in (38) has to be 


replaced by an invariant four-dimensional approach, which will briefly be discussed in the 
next Section. 


2.3 HBT effect for binary stars 

In the previous sections it has been shown that classical electrodynamics and quantum- 
electrodynamical treatment leads to the same expressions for the correlation function of 
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two pointlike sources, see Eq. (14). Therefore, we can just apply any of these expressions 
in order to derive the correlation function for binary stars, given by 


^Binary 9 42 [ ^ "*T [ ^ 

^ ^-star J 


X 2 


1 + cos k 


r 12 ■ d 


AB 


where the regions of integrations are elucidated in Fig. [5j 



(41) 


Figure 5: The both components of the binary system are denoted by 
figures elucidate the region of integration in Eq. (41). 




The integration in (41) yields for the correlation function for binary stars the following 
expression: 


C*Binary B ) — 1 T 


4 L 2 


k 2 d\ B R 2 


T 2 


'k Rd 


AB 


1 + COS 


k a (t) dAB 


(42) 


which is valid for a(t) >2 R ; the extreme case a(t) < 2 R is not relevant for our inves¬ 


tigation. The term in (42) which is not proportional to the cosine-function agrees with 


(38), and describes the standard HBT effect for stars, that means photons from star A 


are correlated with each other and photons from star B are correlated with each other. 


The term in (42) which is proportional to the cosine-function describes the HBT effect 


due to the binary system, that means photons from star A and photons from star B are 
correlated with each other. 

Two configurations are on the scope of our investigation. First, the case of inclination 
i — 0 , where for circular orbits (eccentricity e = 0 ) the apparent distance between both 
stars as seen from the observer is simply given by 


a (t) — A — const. 


(43) 


Second, the case of edge-on binaries (inclination i = 7 t/4 ), where for circular orbits the 
time-dependent distance a (t) between both stars as seen from an observer takes the form 


a (t) 


A cos 



for a (t) > 2 R. 


(44) 
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Here, t is the time, while T is the orbital period, given by 


T = 2vr, 


A 3 


\ G (M a + Mb) 


( 45 ) 


where A is the semi-major axis, and G = 6.67 x 10 


-li 


m 


kgs 2 


is the gravitational constant. 


We will show that these results for the correlation function for stars (38) and for binary 


stars (42) can be reproduced by means of the more sophisticated approach of differential 


HBT-method, which has originally been introduced for analyzing the HBT effect in the 
fireball of heavy-ion collisions mmm- 


3 HBT in heavy-ion-physics and astrophysics 


By means of the HBT effect in astrophysics one may determine the spatial size of a star 
or the distance between the components of a binary system, while the HBT analysis in 
heavy-ion physics allows to get information about the space-time structure and evolution 
of the fireball created in heavy-ion collisions. But both approaches differ significantly. The 
HBT approach used in heavy-ion physics is more involved and more general than in case of 
astrophysics, simply because the hot and dense hadronic fireball expands rapidly in time 
and space on a timescale of about 10~ 22 ... 10 -21 seconds, while a star emits continuously 
thermal radiation and the star’s surface does not alter significantly over a very long period 
of time of about 10 16 seconds. Therefore, a simple two-dimensional integration procedure 


like in (16) is inapplicable in heavy-ion physics, instead a four-dimensional description is 


necessary. Both appraches are compared and the conditions are elucidated about how to 
modify the HBT analysis used in HIC for the case of astrophysics. 


HBT analysis in heavy-ion physics: Since the fireball expands rapidly, in heavy-ion 
physics the second-order correlation function in (16) is generalized into a four-dimensional 
description in momentum-space, see e.g. [32], where it reads: 


C(k A ,k A ) 


P 2 (k A , k B ) 

Pi (k A ) Pi (k B ) ’ 


(46) 


where k A = {oj k A , k A ) and k B = ( oj kB , k B ) are the four-momenta of the two photons and 
P\ and P 2 are the inclusive one-photon and two-photon distribution function, respectively. 
That means, P\ (k A ) is the probability of detecting one photon of momentum k A in de¬ 
tector A or B, Pi (k B ) is the probability of detecting one photon of momentum k B in 
detector A or B, and P 2 (k A , k B ) is the probability of detecting two coincident photons of 
wavenumbers k A and k B in detectors A and B. Let’s assume the photons are emitted at 
two points on the fireball’s surface, X\ = (p, Xi) and x 2 = (t 2 , x 2 ). Then, the one-photon 
and two-photon distribution function is then given by 

Pi (k A ) = J d A x 1 S (xi, k A ) , Pi (k B ) = j d A x 1 S (zi, k B ), (47) 

P 2 {k A ,k B ) = j d A Xid A x 2 S( Xl , k A ) S {x 2 ,k B ) |Ti 2 | 2 , (48) 
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where Ti 2 is the two-photon wavefunction in (34) and its absolute square in (35), and 


the functions S (aq, kjC) and S (x 2 , ks) are the socalled photon emission functions, a term 
originally introduced in heavy-ion physics and which will be considered in the next Section. 

I 11 heavy-ion collisions the HBT effect was originally been used for the same purpose as 
in astrophysics, namely the determination of the system size by measuring the correlations 
of bosons (pion pair correlations and later photon pair correlations) which are emitted 
from the expanding fireball created by the heavy-ion collision [3]. Up to now, the HBT 
analysis remains the only experimental and model-independent approach to determine 
the size af the expanding fireball. Ever since, the effect of photon (or pion correlations) 
has become a fundamental aspect not only in field of quantum optics, but also in the 
theory and experiment of heavy-ion (or nucleon-nucleon) collisions. That means, similar 


second-order correlation functions (46) have also been measured in heavy-ion collisions 


for pion-correlations [H [5] (where T 12 would then be the two-pion wavefunction and the 
four-dimensional integrations would run over the fireball). 

Todays highly sophisticated high energy heavy-ion experiments can even detect the 
ellipsoidal shape of the fireball and also its tilt by means of the HBT effect j6j |33]. At 
present collision energies the angular momentum and rotation of the hot and dense fireball 
is becoming a dominant feature of heavy-ion collisional processes, and the possibility to 
analyze such effects by the so-called differential HBT method was subject of the investi¬ 
gations in mmm- 

In heavy-ion reactions numerous particles are registered, of the order of thousand, 
in every single collision event. The timespan of the reaction is of the order of lOfm/c 
just as the size of the system. So, the particles emitted and observed in one collision 
event are considered as contemporary, and interacting with each other. Consequently an 
event by event measurement registers an M-particle correlation (where M is the bserved 
multiplicity of the event). For the second-order correlation function every pair is selected 
from the M particles and this provides the probability distribution P 2 (Aqq &b) which allows 


to determine the second-order correlation function in (46). 


HIB analysis in astrophysics: In order to exploit the advanced machinery of HBT 
analysis developed in heavy-ion physics, in our subsequent investigation we will start to 
consider the correlator in (46) instead (16). The results of such an approach will finally 


yield similar expressions for the correlator as given by Eq. (38) for stars and by Eq. (42) 


for binaries. I 11 doing so, several aspects have carefully to be accounted for. 

First, as mentioned above, in case of stars and binary stars we face a continuous 
emission of thermal radiation in time. We have to assume that the luminosity from 
the stars is big enough, so that in a correlation time period we have sufficient amount 
of photons emitted in the direction of the observer so that the observable multiplicity, 
M 4. According to the investigations in [34] we can savely assume that this requirement 
is obeyed for many stars located in the neighborhood of the Solar system; e.g. from Sirius 
10 6 photons/cm 2 sec arrive the Earth’s surface. Then, the subsequent analysis can be 
repeated in every such time interval, similarly to subsequent collision events in heavy-ion 
reactions. However, in case of faint and optically not resolvable binaries, the luminosity 
might not be too large, so it is also important that the detection efficiency should be as 
high as possible. In this respect we mention that most modern optical technologies allow 
to generate CCD detectors with a plane of about 1 square — meter, as used for instance 
in the ESA astrometry mission Gaia 
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Second, we have to account for the fact, that the opacity of a star is so extremely high 
that the visible photons do not originate from inner regions of the star but from the thin 
photosphere of the star, that is to say from the two-dimensional surface. Accordingly, the 
HBT approach used in heavy-ion physics has to be modified for the case of astrophysics, 


by simplifying the four-dimensional integrals in (47) - (48) by integrals which run over 
the two-dimensional hyper-surface of stars. 

Third, the wavefunction in (34) is valid both for the geometry of HIC (i.e. R 


and R L) and for the geometry of astrophysics (i.e. d ab N R and R <C 7), and will 
be simplified and designed for the specific case of stars, cf. considerations in Section [2} 
These issues will be the subject of Section [4] and Section [5j 


4 The Emission Function 


Let us consider the photon emission function S (x, k) which is a fundamental ingredient 


in the correlation functions in (46) via Eqs. (47) - (48). The photon emission function, 


frequently called sourcefunction, is the probability distribution of emitting a photon from 
the source point x M = ( ct,x ) with four-momenta k 1 ' 1 = ( k° , k), that means the number of 
photons, AN, emitted in the phase-space element A 3 x A 3 k per unit time At. A Lorentz 
invariant scalar can be obtained by multiplying the photon-energy k° = u>k of the emitted 
photons: 


S(x,k) = k c 


AN 


At A 3 x A 3 k 


(49) 


Accordingly, the total number of photons is given by the summation over the entire space- 
time and momentum-space: 


N = 



(50) 


For the calculation of the photon emission function in (49) one needs the photon four- 
current, defined by 


N»(x) 


/ 




^t/M, 


(51) 


where f{x, k) is the invariant scalar phase-space density distribution (black body distri¬ 
bution) of the emitted photons, 


f {x, k) 


9"f 1 

( 27 t) 3 exp ( k^ ujT) - 1 ’ 


(52) 


where the degeneracy-factor g 7 = 2 accounts for the two polarizations of the photons. This 
single particle distribution depends on the local four-velocity u^^x) of the light-source. In 
order to obtain the total number of photons crossing a space-time hyper-surface E, we 
have to integrate over all momenta and over the entire surface: 


N 


^-k^ffak) 






(53) 
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where dT , M is an infinitesimal three-dimensional surface-element. By not performing the 


momentum integral in (53) we recover the invariant triple differential cross section by the 
Cooper-Frye formula [57] : 

,dN 


dl\ r 

e m = L f(x ’ k) 


(54) 


with E = k°. On account of Eqs. (50) and (53), we get an integral relation between 


photon sourcefunction S (x, k) and photon density distribution f (x, k ) as follows: 


J d‘x S(x,k ) = I ( f(x,k) i J ' dS „. 


(55) 


In the r.h.s. of Eqs. (53) - 


it is assumed that the emission takes place through a 
three-dimensional hyper-surface £ with the normal vector dT,^ = d£ M (x). The hyper¬ 
surface is bordered by the outgoing light-cone, so all emitted particles should cross the 
hyper-surface. We are interested in the hypersurface x° = const.. Then, for time-like 
surface we have = dx dy dz d^, with d^ being a dimensionless unit four-vector. We 
will assume the simplest case for the time-like normal d^ = (1,0, 0,0), but actually we 
do not need this requirement to determine correlation function since the invariant scalar 


Ad 1 dn cancels in the normalized correlator. Then, according to relation (55) and in virtue 


of Eq. (52) the source-function can be parametrized as follows; cf. mm 

2 7s 


S (x, k ) = 


(27t) 3 exp (AT Uf, (x) /T (x)) - 1 


H(t) G(x,y,z) AT <7 M (x) 


(56) 


where % = \Jl — is the Lorentz factor with v s = 


v s 


being the absolute value of 


spatial velocity of the source; let us recall that the chemical potential for photons is zero 
(i.e. fugacity factor equals 1). The function H (t) governs the time-dependence of photon 
emission, which in case of heavy-ion collisions is usually described by a delta-function 
(sudden freeze-out) or an exponentially decreasing Gauss-function (gradually freeze-out), 
e.g. HU, because the fireball emits photons during a small time-intervall around the freeze- 
out time t 0 . Here, in astrophysical systems, the thermal light-source (star, binary star) 
permanently emits photons, that means the number N of emitted photons is proportional 
to time t , i.e.: N ~ t ), hence the function H (t) in (56) remains basically a constant over 
a very long period of time, thus differs considerably from the case of HIC scenarios. The 
function G (x, y, z) is the space-time emission density across the layer of the hypersurface. 
For stars the layer is narrow, i.e. the diameter of the photosphere of a solar-like star is 
much smaller than the diameter of the star, hence the function G (x, y, z) can be simplified 
by a two-dimensional surface and diameter of the layer in z-direction is described by a 
delta-function. Using a Gaussian profile for a star with mean-radius i?, we have 


G(x,y,z ) = 5 (z) exp 


x 2 + y 2 ' 
2 R 2 


(57) 


and the Boltzmann approximation (Jiittner-distribution), we finally arrive at the source- 
function for a two-dimensional surface of a star with radius R: 

S (x, k) = exp (~ k ) H ^ 6 ^ 6XP [^ X 2TR ) ^ ’ ( 58 ) 

where C 7 = 4 7 t 3 . In what follows we will use this expression for the photon emission 
function. The generalization for the case of two light-sources is straightforward. 
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5 The Correlation Function 


The correlation function in Eq. (46) is defined as the inclusive two-photon distribution, 


normalized by dividing with by the product of the inclusive one-photon distributions. 
The correlation function in momentum-space depends on k A and k B being the photon 
four-momenta detected at detector A and B, respectively. We will assume photons are 
emitted at two points, X\ and x 2 on the stars’s surface, which are separated in space, but 
due to the distant source we cannot resolve where the photons are coming from. 

5.1 The wavefunction for stars 


In order to determine the two-photon distribution in (48) one needs to implement the 
two-photon wavefunction in (34), which can be written as follows]^] 


*12 = 


V2 


J[ki A -(x A -xi)+k 2 B -(x B -x 2 )\ _|_ Q i[k 1 B (x B ~xi)+k 2 A (x A -x 2 ) 


(59) 


Recall, the Boson wavefunction in (59) is fully symmetric, irrespectively if we symmetrize 


for xa <->• x B or Xi «->• x 2 , and is valid for heavy-ion collisions and stars. We can consider 
that k = \{kiA + k 2B + fere + k 2 A}- Then the three-momentum vectors become (note 
that these vectors are not independent of each other: k\ A — k\ B = k- 2A — k 2B ), 


ft 


where 


kiA = k - — + — 


ft 


k lB = k - — - - 


i _ , , K <1 

k 2B - k + - - - , 

ko A = k + - + - , 
2 2 


x x -x 2 r 12 . 

ft = k ---= —— and q = k 

1j 


xa~ x B _ ^ d AB 


L L L L 

In case of stars we have d AB r i2 , and the four fc-vectors can be expressed as: 


(60) 

(61) 

(62) 


kiA 

= k- 

ft 

2 ’ 

f^2B 

= k -|— , 

2 

(63) 



ft 


ft 

(64) 

k\ B 

= k- 

2 ’ 

k 2A 

= k+ . 

2 


With these parameters the wavefunction becomes 


*i 2 ar = -j= exp [(ik ■ (x A + x b -Xi- x 2 )] 


x exp 


ft 


-i — ■ (XA-Xs-Xt+X^ 


+ exp 


ft 


+i j ■ (xi-x 2 +Xa-x b ) 


(65) 


Now in order to be able to perform the integrals over the source variables, X\ and x 2) we 
insert the variable ft = k{x\ — x 2 )/L and we describe the remaining observable parameters 


5 V / 7 t \ | 7 ® A 7 7 **T 7 7 *^2 7 7 ^ B '*'2 

°where (ki = k 2 ): k^A — k -—-, kis = k -—-, k 2 A = k -—-, k 2 s = k -—-. 
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in terms of q as d AB = Lq/k , and then we get 


.Star 


l2 = — exp [ik ■ (x A + x B -x i - x 2 )\ 


x I exp 


i k 
2L 


(x 1 - x 2 ) ■ (^y-(xi-x 2 )j + exp + ^(aq - x 2 ) ■ (^Y+(x 1 - x 2 )^j 

( 66 ) 


which, after performing the multiplications in the exponents of the last two terms, becomes 


^Star = - 


1 


12 


V2 


exp 


x 


exp 


ilk - (x A + x B —x i - x 2 ) + jixi - x 2 ) 2 


+ exp 


+M 2 ' 


-i(~- (Xi-Xz) 


(67) 


The structure of this function is the same as of Eq. (34), and then for the absolute value 
square of the wavefunction we obtain: 




Star 

12 


2 1 

= 1 + - ( exp [+iq ■ (x 2 - aq)] + exp [—iq ■ (x 2 - * 1 )] 


= 1 + cos k 


d A B ■ r 12 ' 


( 68 ) 


where we have used the relation q = kd AB /L. The result in (68) agrees with the earlier 


obtained square of the wavefunction in Section 2.2 cf. Eq. (37). Thus we express k in 


terms of the initial state positions k = fc(aq — x 2 )/L, and we perform the integral over 
the source points included in n also. 

5.2 The correlation function 


With the result for the absolute value of wavefunction in (68), we obtain for the one- 
photon distribution: 


P\(k) = j d 4 x S(x,k), 


(69) 


while the two-photon distribution becomes: 


P 2 [k+ ^ - | - |) = J d 4 x 1 cl 4 x 2 S (^c 1 ,k + ^ S (x 2 ,k 


kx 2 q' 

Yl ~ 2 


x 


1 + - (exp [iq ■ ixn - as,)] + exp [-iq ■ (x, - as,)] 


. (70) 


Using Eqs. (69) and (70) with the wavefunction in Eq. (68), together with the definition 

R(k, q) 


of the correlation function, we have: 

C(k,q) = 1 + 


|J d A x S{x, k) | 


2 > 


(71) 


22 































where 


R(k,q) = Re [j {+) (k,q) J ( _)(fc,-g) . 

Here R(k, q) can be calculated via the function 

J(±)(k,q) — I d^x S(x,k(l±x/2L)) exp(iqx) 

= I d 4 x S(x, k{l±x/2L) [cos(ga;) + i sin (qx )], 
and we obtain the R(k, q) function as 


(72) 


(73) 


R(k, q) = j'fx, (.txi cos. (q ■ ( x i - x/j j S j S (x 2 , k - (T . (71 J 


This can easily be verified, by using Eq. (73), forming a double integral over d A Xid i x 2 


from J( + )(fc,g) J(-)(k, — q), yielding to a term exp[— iq(xi — x 2 )\- Then taking the real 
part of the double integral leads to a term cos[g(a;i — x 2 )\ and this recovers Eq. (70). 


5.3 Source with Black Body Jiittner-distribution 


Let us consider the S(xi, k\) S(x 2 , k 2 ) term in Eq. (70). According to Eq. (58), we as¬ 


sume that the single photon distributions, f(x,k), in the source function are Jiittner 
distributions, which depend on the local velocity, u^^x), via the term: 


exp 


—k^ L u^x) 

W) 


(75) 


For a normal star the surface temperature is of order 6000/1, and for optical light in nano 
meter range the exponential is of order 10 2 , so that the Jiittner distribution is a proper 
approximation. In general, the local flow velocity might be different in different locations, 
X\ and x 2) and this fact influences the correlations of the observed momenta. Thus, the 
scalar products in terms of k and hi are: 


-fci «1 p -f!2 U2 _ c -(fc+K/2)«i p -(fc-K/2)«2 _ p-fctll p -fcll2 p-KUl/2 p +K«2/2 


= e 


(76) 


where we used the notation u\ = u(x i) = u^(x i). We assume that for a given detector 
position the normal direction of the emission is approximately the same, so for the two 
sources the term k^ l a^(x) is the same and it cancels in the nominator and denominator. 


Thus, the expression of the J-function in Eq. (73) will be modified to 


J(±){k,q) = / d 4 x S(x,k ) exp 


± 


k x ■ u(x) 
2 LT(x) 


exp (iq ■ x) , 


(77) 


and subsequently we can calculate R in Eq. (72), 


R(k,q) — j d A X\ d i x 2 S(xi, k)S(x 2 , k) exp 


K\ U i 


n 2 u 2 


2 T(x 1 ) 2T(x 2 ) ) 


cos [<7 • (aq - x 2 )} , 
(78) 
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where K\ = k(x±),K 2 = k(x 2 ), n(x) = kx/(2L ), and with the correlation function in 


Eq. (71). The Eqs. (77) and (78) are consistent with the definition (72) of R in terms 


of J(±)- One has to keep in mind that the integrals in Eq. (78) are extended over K.(aq) 
and K.(aq) also. If we have few point like sources the integral becomes a sum over the 
sources, and then the k* values should be taken at the same position as the arguments of 
the velocities, u(xj). 


6 One Source 


6.1 One source at rest 

We will determine the correlation function for one source as given by Eq. m- First 
we consider the invariant scalar k^u^, which can be calculated in the frame where the 
surface-element is at rest. We have then tA = (1, 0, 0, 0) for the four-velocity, hence 

k^u 11 = k°, (79) 

with k° = Ek being the energy of one photon-mode in the rest frame of the star. According 


to (58) and (79), the source function for one source at rest reads 

Ek 


S(x,k) = ^ exp t 


H it) 5 (z) exp f- X ) k 1 ' & fl (x ), (80) 


where T s is the temperature of the source. The denominator in ( ffT| ) is the single photon 
distribution for which we obtain: 


p+OO 


d A x S(x, k) = (ATdyJexp (—7^) I e dx I 

\ 1 X / *7—00 J—O O 


+°° _ y 2 

e 2 Wdy 


D 


= 2 7r — {k^a^) exp ( — — ) , 


Ei 


(81) 


where we have defined D = J dt H (t) (in the normalized correlator (71) this factor cancels 
out), and we have used f dz S (z) = 1 as well as: 


r+00 


■~ ax2 dx = J- 


(82) 


In this simplest case we also assume that the surface direction is a 11 = (1,0, 0,0) where 


the observer is located in the z-direction. The nominator in (71) is determined by means 


of Eq. (77), and we obtain 


J(±)(k, q) = d 4 xe iqx e 


4 , 7 . e -F K °/{2T S ) 


S(x, k ) 


n eu±kP /2 r +00 

(k^)e - 


p+OO 


c. 


D 


e zi&e iqxX dx 

— OO J — OO 

0 


= 2 7T R 2 — (k^dfj) exp 


Ek 


exp 


K 


T 


2 To 


exp 


e 2 B? e iqvV dy 
R 2 


-q 


(83) 


where we used exp(— p 2 x 2 ±qx)dx = (vW p) x exp(g 2 /(4p 2 )) [40]. 
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In the J(k, q)(+) J(k, — q)(~) product the terms exp[±ft°/(2T s )] cancel each other. In¬ 
serting these equations into (71) we get for the correlation function for one source 


C(k,q) = 1 + exp (-.R 2 g 2 ) . 


(84) 


The correlation function for one source at rest does not dependent on k. One might 
want to compare the correlation function in ( [84] ) with Eq. ©• However, one should 
not wonder about the slight difference, because in order to obtain (84) we have used a 


Gaussian source function in (57), while in ( Jl7| ) the Einstein-Hopf model for the thermal 
radiation source has been applied. 



d[m] 


Figure 6: The correlation function for one source in Eq. 

Gaussia 

10 6 m-\ L = 4.14 


(solid line) and the correlation 


function or two spherical Gaussian sources in Eq. (94) (dotted line). The parameter are: 
R = 1.2 • 10 9 m, k = 7.85 


10 


16 


m, x < 


= ar. 


As Fig. [6] shows, the correlation for such a star size yields an extended distribution in 
terms of the detector distance d^B, which corresponds to qAB as qab = k dAB/L. Because 
of the symmetry of the source the correlation function only depends on the absolute value 
of Qabi s ° there is no preferred direction in the plane perpendicular to the source. The 
dependence of the source distance to the correlation function through qAB limits the region 
where the correlation function will be applicable. Systems under investigation are bound 
to the near-zone of the Solar system. 


6.2 One source in motion 

Let us now consider one source which moves in the z-direction with a velocity v z . Then 
we have, u% = 7 S (1, 0, 0, v z ), and the scalar product k ■ u s /T s = k^u^/T s provides an 
additional contribution to the correlation function. However, in the case of a single star 
the velocity and the temperature do not change within the star, so the modifying term 
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in Eq. (78) becomes unity, and we have 


kfi id (E7 k z u 2 ). 


According to (58) and (85), the source function for one moving source 


S (x,k) = H (t) 5 (z) exp 

Ozy 


% (E k - k z v z y 

A 


exp 


x 2 + y 2 ' 
2 R? 


(85) 


k^d^ix). ( 86 ) 


Within the source the velocity u s and temperature T s are assumed to be the same. The 
spatial integrals can be performed in the rest frame of the source, giving the same integral 
result as above (81), because the moving cell-size shrinks, but the apparent density in¬ 
creases, so that the total number of photons in a cell remains the same as it is an invariant 
scalar. For the integral of the one-photon contribution we obtain 


J d 4 xS(x,k) = 2irR 2/ y s (k^a^) 

The two-photon distribution results in 


D 

c. 


exp 




(87) 


J{k,q){±) = J d A xe iq ' x S(x, k) exp 


k ■ Uc 


T- 


D 


= 2nR^ s (k/ l d tl ) —exp 


2 T s 
k ■ u s 


exp 


K ■ Ue 


T- 


2 T„ 


exp 


~q 2 ) ■ ( 88 ) 


When calculating R(k, q), in the J(k, q)( + \J(k, —q)(~) product the terms exp[±K-u s /(2T s )] 
cancel each other. In the formulae the k and k are considered as the wavenumber vectors. 

We then insert these equations into equation (71) and we get for one moving Gaussian 
source 

C{k, q) = 1 + exp (~R 2 q 2 ) ■ (89) 

Again, this result does not depend on k, just as the previous single source at rest in 
Eq. p4|). 



Z 



Figure 7: Two steady sources with a distance between them of 2d in the a:—direction. 
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7 Two Sources 


7.1 Two Sources at rest 


For emission from two steady sources in Heavy-Ion Collisions, two-photon correlations 
were studied in Ref. [55] . Here we use our represented method. We assume that the two 
source system is symmetric, i.e. both positions of the sources are placed symmetrically 
and also their normal vectors, cF, are the same, Fig. [7j If the normal were = (1, 0, 0,0), 
then the invariant scalar k M <3^ = k°] it has been mentioned already that we actually do not 
need this additional requirement to illustrate the correlation function, which would arise 
from an idealized symmetric system. In case of two symmetric sources at rest (y s = 1) 
the source function is (for stars we can savely assume Hi ( t ) = H 2 (£)): 


S(x, k ) 


E s.{x,k) 

s= 1,2 


(k"a„)H(t) E G ‘ {X ' V ' Z) 


s=l,2 


exp 


k ■ u s 
~T 

S 


(90) 


while the function J in the Jiittner approximation is 


J(±)(k,q) 


E ex P 

8 = 1,2 


hi • U s 


exp (iqx s ) / d A x S s (x,k ) exp (iqx), 


(91) 


where x s is the position of the center of the source, and the spatial integrals run separately 
for each of the identical sources, i.e. we assume stars with identical density profiles, but 
with different temperatures, T s . 

In case of steady sources u% = (1, 0, 0, 0), and the spatial integral for one source is the 
same as for a single source. Thus, 


J d^x S ((x, k) 



( 2 vri? 2 ) (kJ'd^) J 2 ^exp 

s=l,2 



(92) 


and 


J(±)(k,q) = ex P 


8 = 1,2 


K 


T 


2 T„ 


(2t tR 2 ) (k^) 


exp (iqx s ) / d A x S s (x, k) exp(iqx) 
•J s 

( R 2 


exp 


q 


D 


x 


E yr ex P ( — jT ) exp 


r 

s=l,2 


hi 


T 


2T b 


exp(i q u x u s ) exp (—iq ■ x s ). (93) 


In the J(k, <?)(+) J(k, — q)(~) product the terms exp[±K°/(2T s )] cancel each other. Both 
J(k, q)(+) and J(k, —q)(~) includes a sum [exp(iq • x s ) +exp (—iq ■ £E S )], and their product 
leads to a factor 2[1 + cos(2qr ■ x s )]. Consequently, if the two sources have the same 
parameter, just different locations, x\ = —x 2 (see Fig.84) then 


C(k, q) = 1 + - exp(— R 2 q 2 ) [1 + cos(2 q ■ £c s )]. 


(94) 
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Like above in case of one source, one may compare the correlation function for binaries in 


(94) with the previous result Eq. (42). The difference between these both results is caused 


by the Gaussian source function in (57) which is used in order to obtain (94), while in 


(42) the Einstein-Hopf model for the thermal radiation source has been applied. 

(Section 9.1 on p.41 ibid.), and in the limit of x s = 0 
If the distance of the two sources is 2d, i.e. 


This result agrees with Ref. 
it returns the single source result, Eq 
X\ = d and X2 = —d, then 2 qx s = 2 q x d, thus the modification appears in the ^-direction 


(84) 


only. In the other directions, q y , the single source result (84) is returned. For distances 
under 4 R the zero-points are artificial because the stars will overlap. 


7.1.1 Sources with varying distances 


Eq. (94) is valid when the two stars are close to each other relative to us. That means, 


they have just past each other in the orthogonal plane. We parametrize the distance 
between the stars, x s with a time component, x s —>■ x s (t). 
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Figure 8: Four diagrams of the correlation function in Eq. (97). The two figures at the top 


show the oscillation of the correlation function over periods 1 a.u. and 25 a.u. respectively. 
The two figures at the bottom show the reduction of the oscillation. The x-axis is to get 
a perspective of the time scales. The parameters are: R = 1.2 x 10 9 m, k = 7.85 x 10 6 m -1 , 
L — 4.14 x 10 16 m, d = 1 m, major axis A is 1 a.u., orbital period 0.7 year. 


Assuming the stars rotate each other in circular orbits, the distance x s (t) on the major 
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axis is given by 


x s (t) =Asin(^t + $), 


(95) 


where A is the maximum distance of the major axis, and T is the total orbital period. $ 
is a phase shift to be adjusted such that when t — 0, the two stars have just past each 
other, x s (0) = 2 R, where R is the radius of the star. Our simple model is not applicable 
when one star is behind the other, shadowing the light from one of the stars. This will 
give us the possibility to have the detector distance fixed, and Eq. (94) will depend on 
time. For two sun-like stars, with A=25a.u., G is Newton’s gravitational constant, the 
period is given by Kepler’s third law: 


T = 2tt, 


A 3 


\ G(M a + M b ) 


= 2.7 x 10 9 s ~ 85 y, 


(96) 


where M A and Mb are the masses of the two stars. Hence, Eq. (94) becomes 


C(q,t) = 1 +-e R2q2 [l + cos(2qx s (t))] 


= 1 + \-e R2q2 [l + cos(2(p4sin(^T T $))]. 

^ / 


(97) 


At d — 1 m there is a rapid oscillation, but there is a rather special behaviour when the 
sign of the gradient of the sine changes. Four figures for two different orbital periods (and 
hence major axis) are plotted in Fig. [8j one for A = 1 a.u., and the second for A = 25 a.u. 
respectively. 

As seen in Fig. [8]there are rapid oscillations, but there is an interesting effect appearing 
during a period of rotation. During the turning points of the sine function the correlation 
function get a reduction in the oscillation. The reduction appears when the stars reaches 
their maximum separation on the axis parallel to the rotation. If the stars are rotating 
in elliptical orbits, the reduction will be different depending if they are at a maximum, 
or minimum separation. Finally, the correlation function over a longer period and over 
different distances is shown in Fig. [9j 


7.1.2 Sources with elliptical orbits 


A binary system with an elliptical orbit, Fig. 10, is characterised by the eccentricity e. 
This will effect the relative velocities they pass each other as seen from our point of 
view, if they are transiting each other when they are closest, or furthest away. In this 
ideal scenario we will lie on the major axis. Then the stars will move faster when they 
are closest to each other. If we imagine that the two stars pass each other when they 
have their furthest and closest distance, to first approximation then they should show two 
slightly different correlation functions after the crossing because of the different velocities. 
We can consider the difference (differential HBT method), 


SC(k,q,t) = i (cos (lq sin (j^i + *,)) - cos Sill ((q + * 2 ))) e " V , 

(98) 

and this function vanishes for e = 0 when they are orbiting in perfect circles; here T\/ 2 
is the orbital period for the furthest and closest distances, and $ 1/2 the two phase shits 
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Figure 9: The correlation Eq. (94) plotted against time t, and detector distance d, 
with the corresponding parameters, distance to system L = 4.14 ■ 10 16 m, wave-vector 
k = 7.85 • 10 6 , phase shift <h = 1.64 ■ ICE 4 , max major axis amplitude A = 25AU, mean- 
radius of stars R = 1.20 • 10 9 m. Second picture is over a turning point. 


respectively and the dependence on the distance with the eccentricity is given by r max = 
A{1 + e) and r min = A{1 — e) [Mj. The period are based on circular motion at these 
distances, so they are only an approximation for a short time. 


X 



Figure 10: Two Stars orbiting around a common center-of-mass in an elliptical orbit. 


Eq. (98) has a very particular shape that can be seen in Fig. Hr and Fig. 11 3. It will 


show its effect in just a couple of days. The closer the system is in a circular orbit the 
longer their correlation functions stays in phase. This can be a method of assisting in 
determination of the orbital period, velocity and eccentricity in case these two stars can’t 
be resolved separately. 


7.2 Two sources in motion 

We study the system the same way as before, but now we use the present method. The 
two sources are moving in opposite directions, so that u s = U\ or u 2 where u'i = % (1, iq), 
u 2 = Ug = ( 7 S ,7<s(—iq)), and u s = 7 s iq, so that U\ = — u 2 , see Fig. 12. Similarly, 
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Figure 11: The differential HBT correlator as given by Eq. (98). The parameter are the 
following: Left pannel: e = 0.01. Right pannel: e = 0.5. Furthermore, orbital period 
T = lyear, stellar mean-radius R = 1.2 • 10 9 m, wave-number k = 7.85 • 10 6 m -1 , distance 


between stars and detectors L = 4.14-10 
A = 1 a.u.. 


16 


m, 


distance of both detectors (Iab = Ira, 


x s = x\ or x 2 where x% = (t s ,x s ), x% = (t s ,—x s ), and X\ = — x 2 . We assume again 
= (1, 0, 0, 0) for the normal of hypersurface, and t\ = t 2 . For binary stars the change 
of the velocity, and position due to the rotation (which is having a frequency of the order 
of 1/day) can be considered as quasi-static, so the two-photon correlation is not effected 
by this change. If we have several sources then the source function is 


S (x, k) = J2 Ss (x, k) = (. k " a^H (t) Y exp 


s= 1,2 

while the J function is 


n 

3 = 1,2 


k ■ u s 

^TT 


(99) 


J{Kq\±) = Y ex P 

s= 1,2 


K • 11 , 


T- 


2 T, 


exp(fga; s ) / d A x S s (x, k) exp(iqx) , (100) 

Js 


where x s is the 4-position of the center of source s, and the spatial integrals run separately 
for each of the identical sources, i.e. we assume stars with identical density profiles, but 
with different velocities, u s and temperatures, T s . 

The integral for one source is very similar as for the case of one single source. Thus, 


/ d A x S (x, k) — / d A x S s (x, k ) 

^ s=1,2^ S 


D 


(k^a^) (2vri? 2 ) — exp 


k 0 ^js 

T x 


exp 


k ■ u s 
T x 


+ exp 


k ■ u s 

T a 


( 101 ) 


This result returns Eq. (92) if u (f = (1,0,0,0). The function J(±)(k,q ) becomes 


J(±)(k,q) = Y ex P 


s=l,2 


K ■ 11 , 


T 


= (k^ d,,) (27 tR 2 ^ exp 


2 T f 
R 2 q 


S J 

2„2 \ 


exp (iqx s ) / d A xS s (x,k) exp (iqx) 
Js 


7 s D 

Y ^^ ex P 


s=l,2 C s 


k • 11, 


exp 


T- 


K ■ U s 


m 


exp (iqx s ) 
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Figure 12: (color online) Two moving sources in the reaction ([x — z]) plane with a 
distance between them of 2 d in the x—direction. The sources are moving in the directions 
indicated by the (red) arrows. 


= (V <v) (2JTA 2 ) exp (- y9 2 ) ^ exp 


k°l s 


k° 7 S 

1 

H- 

- 1 

exp 

+ 2 T s_ 


om\ 


exp (iq"x\ 


x 


exp 


fc ■ u s 


exp 


K ■ U 


2T, 


exp(—fq • x s ) + exp 


fc ■ u* 


exp 


K ■ U 


T- 


2 To 


exp(zq • x s 


( 102 ) 

where the factor exp(fg°x°) can be dropped if the time distribution is simultaneous for 
the two sources, because then x° s = 0. We have again approximated n as a constant, 
k = kQ. This allows us treat k as a constant of the integrals. This returns Eq. (93) if 
u% = (1,0, 0,0). Now we can divide the two-photon correlation with the square of the 
single photon distribution 

. 2 k ■ u a 

Re[J {+) (k,q)J H (k,-q)] 22 

- 12 - = exp (-it q 


cosh ^ ^^rp Us ^ + cosh cos(2 q ■ x s ) 


1/ d 4 x S(x, k)[ 


cosh 


2 k ■ u. 

To 


+ 1 


(103) 


Consequently, if the two sources have the same parameter, just opposite locations with 
respect to the center, and opposite velocities, then the correlation function is 


C(k,q) = l + exp(—iTg) 


cosh ^ 2 k^ u s ^j _|_ cog ] 1 cos (2 q ■ x s ) 


cosh 


2 k ■ u. 

To 


+ 1 


• ( 104 ) 


This expression returns Eq. (94) if — (1,0, 0, 0), and C(k, q) — 2 if q — 0. 


7.3 Two simplest configurations 

We will now consider two configurations 1 and II, where in case I the orbital, [x, y\ plane 
is perpendicular to the direction of the observer, [z ], on the Earth, and II where the Earth 
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Figure 13: The trajectory of the Binary system is in the \x, r/]-plane, the observer is in 
the orthogonal, ^-direction to the plane. The z-axis is the rotation axis. 


is in the plane of rotation, in direction z. We assume furthermore that the two stars of 
the binary have the same mass, temperature and their orbits are identical circular orbits. 
This is a highly simplified configuration compared to the general configuration in Fig. [l] 


7.3.1 I - Orbital plane is orthogonal to the direction of the Earth 


This configuration is shown in Fig. 13 


Let the angular speed of rotation be cu, and k q, k = (0,0, Ax) but the small 
differences in q can be measured. We need to determine the source parameter, u s and 
x s . If the radius of the orbit is A and the angular velocity of the rotation is oj then 
u s = ojA and x s — d — A. The two stars are opposite to each other so that ua = u s , 
ub = -wa, xb = —xa- These vectors are time dependent: So, for anti-clockwise rotation 
x s = x s (cos (cut), sin (cut), 0 ), and u s = u s (— sin(o;t), cos (cut), 0 ). 

In case I, the vector k is orthogonal to both ua and ub all the time, so we can have 
a z directed k only, and so (k ■ u s ) vanishes. The products (q ■ x s ) and (k ■ u s ) for z 
directed q vanish so 

C{k z , q z ) = 1 + exp (~R 2 q 2 z ) ■ (105) 

Usually R A, therefore C(k z ,q z ) ~ 2 for astronomical configurations, and the depen¬ 
dence on the much larger A leads to stronger variation at smaller q values: 


C(^) = l+iexp(-fl^) 1 + cos (2q x A cos(wt)) cosh 


C(k z , q y ) = l+-exp (-R q y ) 


1 + cos (2q y A sin(o;t)) cosh 


(—n x ujAsm(ujt) 
V T s ~ 

' Kyld A cos (ut) 


. (106) 


Here we assumed that the orbital velocity of the binary star components is non-relativistic, 
so we can neglect the relativistic 7 factor. k x and n y are respectively the mean angular 
size in x, and y direction. For a symmetric source they are both n x = n y = k = kQ. 

As R <C A the leading exponential term tends to 1 for small and intermediate q x or 
q y values. At the same time at q x = q y = 0 both the cos and cosh terms tend to 1 , so the 
value of the correlation function is 2. The cos terms drop to zero when q x = 7 t/[2 A cos(cut)] 
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Figure 14: The trajectory of the Binary system is in the \x, z]-plane, the observer is in 
the plane in the ^-direction. 



> 


Figure 15: The trajectory of the Binary system is in the [y, z]-plane, the observer is in 
the plane in the ^-direction. 


or q y = it/\2Asm(ujt)\. These terms may drop much faster than the leading exponential 
term as A^> R, and show a rather special t depencence due to the orbital rotation. 

The arguments in the cosh terms are of the order ps 10 -8 for a binary configuration. Same 
with the two stars orbiting with a distance 1 a.u.. This is no surprise because the extra 
terms are coming from the velocity dependence of the emission fun ction through Jiittner 
distribution. The velocity is given by cuA ps 400 — . Thus Eq. (106) effectively goes 


back to the result presented earlier when the velocities of the sources where not explicitly 
accounted for. 


Eq. (106) goes like the earlier result Eq. (97), as a consequence of the velocity- 


dependent part of the correlation function going to unity. 


7.3.2 II - The direction of the Earth falls in the orbital plane 

This configuration can in principle occur in two ways, the z-axis is spanning the distance 
between the center of the Binary and the observer on the Earth, and the plane of rotation 


of the binary can either be the [x, 2 :]-plane or the [y, z]-plane, see Fig. 14 and Fig. 15 
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We can freely choose the direction of the x and y axes, so we chose the configuration 
so that the rotation plane of the Binary is the [x, z]-plane. Consequently, with this choice 
neither the distance between the two starts of the binary nor their velocity will have any 
^/-component, see Fig. [T4| Thus in case II the relation between u s and x s are the same, 
but now the anti-clockwise rotation is in the [x,z] plane, x s = x s (cos(aA), 0 , sin(a+)) and 
u s — u s (— smut, 0 , cos(u)t)). 

The decomposition of the correlation function, Eq. (104) now becomes: 

cos(2q x A cos (cot)) 


C' (^20 Q.x) 

= l + exp (~R 2 q 2 x 

C {kzi Qy) 

= l + exp (~R 2 q 2 y 

C(k z , q z ) 

= l + exp {—R 2 q 2 z 


cosh ( 


2k z cjA cos(ujt) 


J + cosh ( 
cosh ( 


—K x uAsin(cdt) 

T s 


2k z ujA cos (ujt) 
Ts 


+ i 


cosh ( 


2k z ujAcos(iJt) 


+ 


cosh cos(2q z Asm(ujt)) 


cosh ( 


2k z u;Acos(ujt) 
T s 


+ i 


(107) 


Just as in the previous case, we have assumed that the orbital velocity of the binary star 
components is non-relativistic, so we can neglect the relativistic 7 factor. Also as we have 
seen the term cosh~ 1, and the argument in cosh (k z u z ) is of order « 10“ 4 . This is 
also 1 for an ordinary binary system. Even for the fastest binary observed with a period 
of 5 minutes the contribution is only of « 1.004 |42j . 


8 Summary and Outlook 

Our investigation starts with a reconsideration of the classical HBT effect, both in terms 
of classical electrodynamics and quantum electrodynamics. Within a simple model for 
the thermal radiation held we have determined the correlation function for stars at rest 
in Eq. ( 0 , and we have explained how one may deduce the diameter of stars from this 
correlator. Afterwards, the correlation function for binary systems has been determined 
in Eq. ( J42| ) for comparsion with a modified HBT approach used in heavy-ion collisions. 
Another aspect which has carefully to be treated concerns the two-photon wavefunctions, 
which differ in case of stars and heavy-ion collisions (HIC), because of the different geo¬ 
metrical limits of detector distances and radial diamater of the thermal sources, i.e. either 
stars or hadronic fireballs. 

Afterwards, inspired by recent developments with HBT analyses in heavy-ion collisions 
taking the relative velocity of the sources into account, we have tried to modify the HBT 
approach used in heavy-ion collisions to the case of binary star systems. Especially, we we 
have considered the differential HBT approach as it is in use for analyzing the hadronic 
fireballs created in heavy-ion collision experiments. 

Since the fireball in heavy-ion collisions expands rapidly in space and time, the cor¬ 
relation function in HIC is determined by covariant integrals over the four-dimensional 
space-time. On the other side, stars emit permanently thermal radiation and ther stellar 
radii remain constant. Furthermore, the opacity of stars is usually so extremely high, 
that the photons only originate from the two-dimensional surface of the stars nstead from 
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the entire three-dimensional volume of the stars. Accordingly, the HBT approach used in 
HIC has to be modified for the case of stars by simplihying the four-dimensional integrals 
itno two-dimensional integrals. 

A further key-point in order to determine the correlator regards the emission function 
which determines the number of photons emitted into a phase-space element. Further¬ 
more, it has been shown that the Jiittner-distribution, frequently used in HIC, is also 
a good approximation for stars. Then the new approach has been applied in order to 
determine the correlator for stars at rest in Eq. (84) and for stars in motion in Eq. (89), 
While an exact agreement between the classical approach and the modihed HIC approach 
cannot be expected because of the different descriptions for the thermal radiation held, 
we have shown that the result in (84) corresponds to the previuos result in Eq. 

Then we have determined the correlation function for binary stars at rest to each other 
in Eq. (94) and for the case of binary stars in motion in Eq. (104). Like in case of one 
star, we have shown that the result in Eq. (94) obtained from the modihed HIC approach 
corresponds to the expression in Eq. (42) which was obtained from the classical HBT 
approach. 

The hrst results of our approach have shown that the dependence on velocities do not 
influence the correlation function for normal binary stars, as one expected. But the rather 
specihc oscillation dependence of the correlation function during a period of rotation is a 
characteristic sign of binary systems, and can help to determine some parameters of the 
orbits when the stars are not optically resolved. This has been elucidated by the difference 
of the two correlators of a binary system (differential HBT approach) in Eq. (98). 

It has been outlined that due to the restriction of real experiments being bounded on 
Earth limits the possibility of fully exploiting the correlation function compared to heavy- 
ion collisions. Our model assumes a thermally equilibrium symmetric sources, while more 
exotic objects could in principle also be investigated within the approach presented. 
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